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The analytical solutions for the third-harmonic generation (THG) on infinite chains in both 
^ ' Su-ShriefFer-Heeger (SSH) and Takayama-Lin-Liu-Maki (TLM) models of trans-polyacetylene are 

^ ' obtained through the scheme of dipole-dipole (DD) correlation. They are not equivalent to the 

results obtained through static current-current (Jo Jo) correlation or under polarization operator P. 
The van Hove singularity disappears exactly in the analytical forms, showing that the experimen- 
tally observed two-photon absorption peak (TPA) in THG may not be directly explained by the 
single electron models. 
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I. INTRODUCTION 

^ ' The nonlinear optical (NLO) properties on conjugated polymers have always been received wide interest for study. 

\ , Prnl Experimentally observed nonlinear phenomenon in conducting polymers, such as photoinduced absorption, 
I bleaching, photoluminescence, ultrafast optical process and exciton behaviors, [||-^ have already encouraged the 

theorists for the possible explanations. |10|-p2[ 

In the early 90s, the large third order susceptibility (x ) of trans-polyacetylene (PA) have been observed exper- 
imentally ||l^,H and have been received wide explanations theoretically. p"5|"|2^ For a long time, the origin of the 
experimentally observed two-photon absorption peak (TPA) |0,n4| receives persistent discussions in theory. |13j24] 
f^ ' Based on the static current-current correlation {JqJo) P^ill^ and followed by Keldysh Green function methods [ p5[ , 
^^ TPA could be interpreted through the analytical forms under both SSH p^ and TLM models p^ for infinite chains. 
^^ Besides the unnecessary application of nonequilibrium diagram techniques in a stable optical response problem, |l9] 
however, the various numerical approaches from dip ole formalism(E • r), such as the results by Yu and Su [|l8yi9| from 
Butcher-Bloembergen-Shen dipole formula, [pSH^OJ the results by Wu and Sun ||2y,gl|] from Genkins-Mednis approach 
"c^ lO'iH ^^"^ *^^ results by Shuai and Bredas from sum-over-state(SOS) Orr-Ward formalism, [Q have shown that the 

^ , TPA is no longer existing if the damping in the energy is considered. ||lq-[24[| Therefore, the two-photon cusp under 
I ' the current schemes is considered to be a van Hove singularity |3J] caused by the singular density of states (DOS) 
''O ' on the Fermi surface in one-dimensional (Id) systems. The fact of the nonresonant property in two-photon cusp is 
^ ■ also noticed by Wu in his current formulations. [ |l5| However, besides the difficulty in explaining this TPA under 
f_\ ] current formulas, there exists another obvious difficulty in theory — apparent zero frequency divergence (ZFD) in 
IL" ■ the definition of the current formalism, pq] This problem actually was realized by Dakhnovskii and Pronin in their 
. ^ , computations based on SSH model. |T^ Unfortunately, the results of x'"^^ obtained from SSH model returns TLM 
J^ ' model if the process of linearization procedures is made under Jq Jo scheme. p6| 

J^ ] Because of the apparent divergence difficulties in current formulas p8| , practically, intense studies based on dipole- 

dipole (DD) correlation have been carried out to obtain the physical pictures of TPA in x^^' of trans-polyacetylene. 
On one hand, the weakly correlated and single electron theories are applied in the numerical dipole computations, thus, 
the finite size effect pi| and the lifetime assumption p^] are seemingly good in generating a peak at exact position 
of TPA; On the other hand, the strongly correlated electron theories have also been applied to short chains, Soos 
and Ramesesha obtained the TPA based on Pariser-Parr-Pople (PPP) models, but it is shifted too low in frequency 



and has too small an intensity p7| ] relative to the experimental peak. |13| Recent work by Guo and Mazumdar 
shows that this TPA could be interpreted as a three-photon resonance based on extended Hubbard Hamiltonian. pj] 
Whether or not the weakly correlated theory or even the single electron models such as SSH and TLM models are 
still suitable in describing TPA still remains as a question. Meanwhile, in obtaining the qualitatively correct result, it 
was suggested by Etemad and Soos that the attention should be made to the x'^' frequency dependences rather than 



magnitudes. |35| It is certainly reasonable in physical point of view since the experiments are more sensitive to detect 
frequency dependences. Thus, some typical approximations are made in the dipole formula. Substitution of the dipole 
operator D by the polarization operator P, namely, the dipole approximation ||2^,Q has been extensively applied in 
the nonlinear optical formula, however, x computations have shown both quantitatively and qualitatively different 
results between P and dipole operator D. [p^ Whether the problem caused by the dipole approximation still existing 
in x*-'^^ is still under our concern. Furthermore, the static dipole formulas |0-p4| and static current formulas [15 28 
are conventionally considered to be equivalent. pq , |2^ -|30| However, as we pointed out in a recent work, p^ it has 
ZFD problem in the static current formalism when the gauge phase factor is ignored, the static formalism under two 
gauges should not return the same results under the exact same basis of unperturbed wave functions. [^ Whether 
this noncquivalcnce between Jq Jo and DD correlations existing in the specific case of x of trans-polyacetylene is 
another interest topic. Although the qualitative features have already been reported in the numerical solutions based 
on DD correlation, ||l^-p3 the exact analytical solutions from the dipole formulas of x'''^^ have not been obtained 
previously. The analytical results, if obtained, will be very helpful for us to directly check all above concerns and give 
more information. 

Fortunately, both SSH and TLM models are single electron models and can be exactly solvable for the nonlinear 
susceptibilities under perturbated schemes, although the efforts will be no doubt paid to cumbersome computations. 
After choosing the long wavelength approximation which are conventionally used previously, [|l5|- p4| , p7[ we could 
directly compared those analytical results between two gauges for both models. The paper is organized as follows, we 
will first discuss the SSH Hamiltonian in Sec.H.A. The analytical form of x*-^-* by the dipole-dipole (DD) correlation is 
obtained and compared with the analytical form by Jo Jo correlation for SSH models in Sec.H.B. In Sec. II. C, we will 
obtain the results under the polarization operator P, the qualitative and quantitative difference between the results 
under D and P are outlined. In Sec. Ill, we will obtain the exact analytical forms under DD correlation for TLM 
models and directly compared with those under Jo Jo correlation. The discussions on the possible implications of this 
work are made in Sec. IV and a brief conclusion is made in Sec.V. 



II. THIRD ORDER SUSCEPTIBILITIES FOR SSH MODELS UNDER DIPOLE FORMULA 
A. SSH Hamiltonian in real and momentum spaces 

Based on periodic tight-binding-approximations. The SSH Hamiltonian P6| is given by: 
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(2.1) 



where io is the transfer integral between the nearest-neighbor sites, A is the gap parameter and CJ ^{Ci^s) cre- 
ates (annihilates) an TT electron at site I with spin s. In continuum limitation, above SSH model will give the TLM 
model. 12^ For the SSH model, each site is occupied by one electron. 

Under the DD correlation, the interaction Hamiltonian is expressed by -ffE-r = — eE • r = — D • E, To provide the 
periodic property and to avoid the ill-definition of position operator f , we should express the position operator f 
under the Bloch states: \n, k >= u„.k(r)e' '"", where n and k are band index and crystal momentum correspondingly, 
Wn,k(r) is the periodic function under the translation of lattice vector. [p8[ Thus, we obtain: 



and 



Tnk.n'k' = «(5„,„'Vk5(k - k') + f7„^„'(k)^(k - k'), 
^n,«'(k) = - / <k(r)VkMn',k(r)rfr, 



(2.2) 
(2.3) 



where v is unit cell volume. 
We change Hamiltonian Eq.(p 



into the momentum space by applying the following consecutive transformations: 
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where 



e{k) = yJ[2tQC0s{ka)f + [Asin{ka)f , (2.7) 

and aj,'^^ (i) and a^"^ (i) are the excitations of electrons in the conduction band and the valence band with momentum 
k and spin s, Ri^ and Ri^ are odd and even position defined by 

Ri=la+{-iyu, (2.8) 

If we choose the spinor description il^l g{t)= ('^fc 5(^)7 'H^sW) ^^^ ^PPly the long wavelength approximation, pq ] 
SSH Hamiltonian including E • r in momentum space is described by: 

HssH{k,t)=Ho+HE.r, (2.9) 

where 

Ho= Y. <k)i;l^{t)a3iJkAt) (2-10) 

and 

iJE.r = -D ■ Eqc^^K (2.11) 

By Eq. ( ^.2[) , the dipole operator Z) could be obtained as follows: [Q 

^ = e ^ (/3(fc)^L^2V'M+^4:ViLV'M). (2-12) 
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where 



is the coefficient related to t he interband transition between the conduction and valence bands in a unit cell 2a and 
the second term in Eq.( 2.12|) is related to the intraband transition, ||3^ e is the electric charge and a are the Pauli 



matrixes. We have already omitted the relative distortion ri{= 2u/a) in this computation since the contribution of 77 
is quite small according to the linear case. [p7|] 

B. Analytical form of x through DD correlation 

Within the semiclassical theory of radiation, |^-Q the electrical field is treated classically and the third order 
susceptibility x''^^ is described by: eF 
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(2.14) 



where 51 = — N cui and T is the time-ordering operator with D is a dipole operator. 

The third- harmonic generation (THG) is defined by setting uji — 102 = ^3- FoUowing simiHar procedures as we have 
done in the hnear computations, p^ by Eq.(2.12) and Eq.(2.14), we obtained the following expression for Xs§h i'-^T^)= 
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where /3(fc) is defined in Eq.(2.12) and hq is the number of chains in unit cross area and the polymer chains are 
assumed to be oriented, and Green Function G{k,uj) defined as following: p^ 



G{k,Lu) 



9 2 i ~ ' 



(2.16) 



with tok = e{k)/h and e — 0+. 

After tedious derivation, we obtained the following analytical expressions for the third-harmonic generation for the 
SSH models rssni^)- 
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The above expression Eq.( 2.17 ) are almost identical to Eq.(A13) in Agrawal et a/.'s work |32| except one more term 
(last term in (A13)), and with the following substitution is applied in Eq.(Al3): 
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Eq.( p.l7 ) is also quite similar to Eq.(ll) in Wu and Sun's result [Eoj based on Genkins-Mednis approach pit] with the 
above substitution (2T8) made and if the symmetric procedures on the negativ e fre quency term should be applied. 



For infinite chains, we could separate two full derivative terms of k from Eq.( 2.17 ) as follows: 
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If we consider following two facts in actual optical process: (i). The velocity on the Fermi surface Vk£{k)\k=±r^ — 0. 
(ii). The life time of the states, thus, we should include the damping in the process, that is, ujk—* ^k + *efe and e^ 7^ 0. 
The above full derivative term It after integral of k then is for the infinite chains. Eq.(2.17) can be simplified as 
following: 
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?iijjfe/A, z = huj/ {2A) and S = A/{2to). F,E,I1 are first, second and third kind 
9| There is exactly no two photon cusp in Eq.( 2.2C| ). This result is consistent with 



, (3) 8 e^no i2toaf 

of complete elliptical integrals 

Yu and Su, Shuai and Bredas' results. And it is quite different from the expression obtained from JqJq correlation 

which shows the following form for Xsshjj ^^ ^^-^ models: |16| 

X^f G^.j. = B z-«{[5 - 8z^(l + 6^) + 2Qz^S^]g(z) - 8[1 - 4^2(1 + ^2^ ^ iQ:,is^]g(^2z) 

+ [3 - 24z2(l + S^) + 18Sz*5^]g{Sz) - 8S^z^E{^, y/T~P)}, (2.22) 



where B = 5xo /(1024z^), same as the constant as defined in Eq.(2.20) and g(z) is defined by Eq.(2.21). We should 
point out that the elliptical form of Eq.(2.22) [Q is the exact the same as the integral form of Eq.(ll) in Wu and 
Sun's result, [Q which was derived from the incomplete Genkins-Mednis approach. |gl|,0 The disappearance of ZED 
under static current schemes is somewhat puzzling or simply should be understood as coincidence for the specific SSH 
model. As we pointed out in a recent work, pq the ignorance of gauge phase factor in the static current formula will 
cause ZED in the susceptibilities and non-equivalent result from the dipole formula. The derivation in the previous 
current work ||l6| does not consider the gauge phase factor's effect, it will cause ZED even for linear case if following 
a straightforward computations. p7| 

Obviously, we observe the qualitatively difference for results under DD (Eq. p.20 )) and JqJq (Eq-(2.22)) correlation 
especially for z — 1/2 and z ~ 1. The inequivalence could be understood by the ignore of the gauge phase factor in 
optical response theory. | p7[ | 

L et X -^ a; + ie in the integral of Eq. (|2.2C| ), the comparison graph of absolute values between Eq.( ^.2C|) and 
Eq.( 2.22 ) for the SSH models are plotted in Eig.l with the following parameters chosen for trans-polyacetylene: |l§H2l|] 
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FIG. 1. Computed DD values (solid line) vs. JoJo values (dashed line) of \Xssh{'^)\ with z = hiu/(2A) and e = 0.03. 

From the graph, we find another very obvious peak at z = 1 at a ratio of 1/10 of the peak at z = 1/3, which is also 
different comparing with the previous theoretical computations. |15|-p4| However, this peak has not been reported by 
the experiments because it is out of the scanning range of photon energy. |13| , |l4| 



C. Results under polarization operator P 



In a recent work, we discussed the different results caused by polarization operator P in linear response p7\ . The 
polarization operator P is extensively applied in the NLO theory. |2^,£3l Based on SSH model, we can do a comparison 
of XssH between D and P. 

The polarization operator P is defined in real space as: 



P^J^E^C/G. 



(2.23) 



If the chain region I is going from 1 to N, we obtain the unit polarization P""'*(fc) from the total polarization 
ptotai^j^^ in the momentum space: 



pumtfj^^ 



hm TT = ^ 2^ V'fe,,ff2V'fe,. 
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Substituting D with P in Eg. ( 2. 14 ), we obtain the XssHp ^'^'^ SSH models as following for infinite chains: 

A56^ n dx f 1 9 1 
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The above expressions lacks TPA. However, the magnitude is too small compared with results under D. {\Xssh I 
from P is about 10"'' of \xssh\ from D ii S = 0.18 for trans-polyacetylene.) Another peak at z = 1/(35) « 1.85 will 
shown up through P • E, this peak corresponding to the transition from the bottom of valence band to the top of 
conduction band. This peak in infinite chains seems not to agree with 'unklapp enhancement' in the solid states, pj] 
While in DD computations, this peak at z = 1/(3(5) will not be very obvious. The comparison on the absolute value 
between Xssh ^^^ Xssh ^^11 be shown in Fig. 2. From the graph, the different shapes of curves are quite obvious. 
Besides the possible ill-definition of polarization operator P, |^,^,Q both the qualitative and quantitative features 
obtained from the P are quite different even for nonlinear optical response theory. Thus it will not be quite suitable 
in using P to obtain the equivalent results as D. 




FIG. 2. Computed \xssHp{'^)\ under polarization operator P (solid line) vs. computed IxsFh ('^)l (dashed line) under dipole 
operator D with z = huj/{2A), e — 0.03 and a magnification of lO'* in \xssHp{'^)\- 

III. THIRD ORDER SUSCEPTIBILITY FOR TLM MODEL UNDER DIPOLE FORMULA 



The results under TLM model for the DD correlation can be obtained by setting (5^0 and (2toa) — > fivp in 
Eq.(2.20). p^ Since we have the following properties for the complete elliptical integrals: [B9| 



f limS(|,VT^^)-l, 
limj2i^(^, Vl-(52) = 0, 
lim n(-, n™, Vl - 5"^) = /(^^z), 
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We will obtain the following expressions for THG under TLM models for DD correlation, 
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where Xo defined in Eq.( 2.2C ). The Xt?.m i^ pfotted in Fig. 3. And no singular property such as van Hove singular 
shows in x^lm under DD correlation. 

In Jo'^ correlation, the result of Xr^Mjj obtained by Wu as follows: |l^ 
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{(5 - 8z2)/(z) - 8(1 - 4z2)/(2z) + 3(1 - Sz^)!{2.z)} 
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The comparison between our result(Z?Z?) and Wu's result {Jo Jo) [|5| on absolute value of Xt^m i^ plotted in Fig. 4. 
If the zero frequency limitation is made, let z — > 0, we obtain the different values for the zero frequency limit for both 
DD and Jq Jo correlations: 
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XTLMiii^) ^ nXo ~ 5.0 X 10 ^^esu. JqJq correlation 



The nonequivalent results between DD and JqJq correlation are quite obvious even for the static limit, 
disappearance of ZFD could be understood as coincidence. |M 
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FIG. 3. The real part (solid line) and the imaginary part (dashed line) of XtlmI'-^) with z = hw/{2/S). 
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IV. DISCUSSIONS 



The analytical forms Eq.( 2.20 ) and Eq.( |3.3| ) show that there will be exactly no two-photon p eak o r even no van 
Hove singularity under DD correlation. They are qualitatively different from the expressions Eq.(2.22) and Eq.(3.4) 
derived from Jq Jo correlation. The results with no TPA peak under single electron models are certainly reasonable 
in physical pictures, which is also agrees with the previous arguments on this problem. [l7[-|24|] The nonequivalence 
between two gauges also provided some evidence that the gauge phase factor in the optical response theory can not be 
ignored, and the static current expressions |15|JlqP^ can not be directly used and return the exact same results as DD 
correlation. On one hand, DD correlation should be much more suitable to study the nonlinear susceptibilities than 
J J correlation. | |37[ | On the other hand, based on the exact expression of dipole operator D, we find that polarization 
operator P can not be directly applied in periodic systems where the electrons are certainly not bounded. And P 
will returns the qualitatively different curves from the rigorous results. From the results we obtained in this paper, it 
seems that the direct application of single electron models may not be suitable in explaining for TPA. It also implies 
that the size effect [Ell, lifetime considerations |l^-p3], exciton effects ftO,E3, quantum fluctuations or disorders 
should be included to explain this non-trivial TPA. Both SSH and TLM models still can be served as the good basis 
to include all those interactions. 



V. CONCLUSIONS 

The analytical solutions to the THG for both SSH and TLM models obtained under the DD correlation com- 
putations, showing nonexistence of the two-photon peak exactly. The nonequivalent results between DD and Jo Jo 
correlations provide an evidence for the possible explanations through gauge phase factor. |3^ In actual practical 
computations, DD will be much more suitable for obtaining reasonable results than Jo Jo correlation, however, the 
replacement of D by the polarization operator P will result in qualitatively different results especially in periodic 
systems. 
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